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Abstract 

We propose an extension of tri-bimaximal mixing to include a non-zero reactor 
angle #13 while maintaining the tri-bimaximal predictions for the atmospheric 
angle #23 = 45° and solar angle #12 = 35°. We show how such tri-bimaximal- 
reactor mixing can arise at leading order from the (type I) see-saw mechanism 
with partially constrained sequential dominance. Partially constrained sequential 
dominance can be realized in GUT models with a non-Abelian discrete family 
symmetry, such as A±, spontaneously broken by flavons with a particular vacuum 
alignment. 
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Over the last decade neutrino physics has undergone a revolution with the mea- 
surement of neutrino mass and lepton mixing from a variety of solar, atmospheric and 
terrestrial neutrino oscillation experiments |1J. Lepton mixing is described by the 3x3 
matrix [2] 
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The Particle Data Group (PDG) parameterization of the lepton mixing matrix (see e.g. 

C12C13 S12C13 si 3 e~ l5 

U— I — C23S12 — Si3S23Ci2e 1<5 C23C12 — si 3 s 23 si 2 e lS S23 c i3 | P, (2) 

S 23 Si 2 — Si 3 C 23 Ci 2 e lS — S23C12 — Si3C23Si2e 1<5 C23C13 

where s\ 3 = sin0i 3 , c\ 3 = cos 9i 3 with 9\ 3 being the reactor angle, S\ 2 = sin^i 2 , 
C12 = cos6'i2 with 9 12 being the solar angle, S23 = sin 9 23 , c 23 = cos ^23 with 6*23 be- 
ing the atmospheric angle, 5 is the (Dirac) CP violating phase which is in principle 

. ai .012, 

measurable in neutrino oscillation experiments, and P = diag(e' 2 , e 1 2 , 0) contains 
additional (Majorana) CP violating phases 0:1,0:2. 

Current global fits for the solar and atmospheric angles typically allow the following 
la ranges (in degrees) [1], 

012 = 34.5° ±1.4°, 023 = 43.1°^:. (3) 

However the situation concerning the reactor angle is more subtle and requires some brief 
discussion. In global fits it has been noted that there is a mild tension between solar 
data from SNO and KamLAND, and a further weak tension involving atmospheric data 
from SuperKamiokande, where both tensions may be resolved by allowing a non-zero 
value for # 13 , corresponding to the following value with la error [3], 

sin 2 0i3 = 0.016 ±0.010, (4) 

corresponding to a 1.2<r indication of a non-zero value for $13. The SuperKamiokande 
collaboration is presently in the process of analyzing its SK-I, SK-II and SK-III atmo- 
spheric data including the effect of #13, so the situation regarding the atmospheric data 
should be clarified in the near future. In the meantime, MINOS has announced its first 
results on electron appearance [6] and observes an excess of events at the level of 1.5a. 
This may be also interpreted as a hint for #13, and this has been combined by Fogli "at 
face value" with the global fit above to give roughly [5], 

sin 2 0i3 = 0.02 ±0.01. (5) 
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Thus, apparently following the adage "many a little makes a mickle", one is led to a 
2cr indication for a non-zero value of #13. This corresponds to a value for i3 in the la 
range (in degrees), 

013 = 8° ± 2°. (6) 

In any case it is certainly theoretically plausible that #13 could take a value in the above 
range [7], so it is interesting to consider this possibility, and we emphasize this more 
general motivation. 

It is well known that the solar and atmospheric data are consistent with so-called 
tri-bimaximal (TB) mixing [8], 
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corresponding to the mixing angles, 

0i 2 = 35.26°, 2 3 = 45 O , 0i3 = O°. (8) 

The ansatz of TB mixing matrix is interesting due to its symmetry properties which seem 
to call for a possibly discrete non-Abelian family symmetry in nature [9] . There has been 
a considerable amount of theoretical work in this direction [101 [TTJ [121 [T3J, [HJ [15]. The 
presence of a non-zero reactor angle as in EqJU] would be clearly inconsistent with the TB 
prediction for the zero reactor angle in EqJH] and so the TB ansatz would be excluded, 
even though the predictions for the solar and atmospheric angles remain acceptable. 

In this paper we shall explore the possibility of extending the TB mixing matrix to 
allow for a non-zero reactor angle 0i3, while at the same time preserving the predictions 
for the tri-maximal solar angle and the maximal atmospheric angle given by EqJHl namely 
012 = 35.26° and 023 = 45°. In order to maintain these predictions requires, 

1. (9) 

To leading order in U e ^ the conditions in Eqj9] correspond approximately to, 

\U e2 \ 2 ^l/3, ~ 1/2. (10) 

We refer to the above proposal as as tri-bimaximal-reactor (TBR) mixing, to emphasize 
that tri-maximal solar mixing and maximal atmospheric mixing are both preserved while 



\Ue2\ 2 


1 




\Uel\ 2 


~ 2' 





1 Note that different versions of the TB mixing matrix appear in the literature with the minus signs 
appearing in different places corresponding to differing choices of charged lepton and Majorana phases. 
We prefer the convention shown which emerges from the PDG parametrization when the angles are set 
equal to those shown in EqlH] 
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non-zero reactor mixing is introduced. To leading order in {7 e3 the TBR mixing matrix 
following from Eq{TU]is then, 
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if* 1 77 \ 
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P. (11) 



V 7^ 1_vz ^ -73^ + 73^ 71/ 

The TBR proposal in EqsJHl EH and [TT] should not be confused with the tri-maximal 
proposal [16] that the second column of the mixing matrix should consist of a column 
with all elements equal to 1/ v^, corresponding to the second neutrino mass eigenstate 
being a democratic combination of flavour eigenstates, i.e. \U[ 2 \ 2 = 1/3 for I = e,fi,r. 
Tri-maximal mixing predicts a solar mixing angle approximately given by 35°, to leading 
order in f/ e3 . However, for |i7 e3 | = 0.14 (corresponding to \U e %\ 2 — 0.02) there are 
significant deviations which tend to increase the prediction for 9u beyond the la range 
quoted in Eqj3] [TTl [18] . In addition tri-maximal mixing also predicts deviations from 
maximal atmospheric mixing which depend on the combination ?7 e3 cos<5 [TTl fT8]. For 
example, according to Fig. 2 of [IB], for \U e ^\ = 0.15, tri-maximal mixing would constrain 
the phase 5 to lie in one of the two approximate ranges 5 = 7i/3—2tt/3ot5 = 37r/2— 97r/5 
in order to yield an atmospheric angle within its la experimental range. 

In contrast, TBR mixing defined by Eqj9] predicts both tri-maximal solar mixing 
9 12 = 35° and maximal atmospheric mixing 6*23 = 45°, for all values of 6*13. However, in 
practice, realistic models only predict approximate TBR mixing, so the leading order 
approximations to TBR mixing in Eqs fTU|[TTl should provide an adequate approximation. 
In addition, as discussed later, model dependent corrections to the TBR predictions for 
the solar and atmospheric angles are expected at some level. Nevertheless, since tri- 
maximal mixing predicts that tan 2#23 is inversely proportional to £/e 3 cos(5 [TTj [TS] , 
while exact TBR mixing predicts maximal atmospheric mixing, these two proposals 
may be distinguished by accurate determinations of 6*23, #13 and cos 5 at future high 
precision neutrino facilities. 

In general, note that TBR mixing differs from all variants based on preserving a 
particular row or column of the TB mixing matrix [18], since the approximate TBR 
requirements \U e 2\ 2 ~ 1/3 and {U^l 2 ~ 1/2, together with U e3 being non-zero, results 
in every row and column of the TBR mixing matrix being different from that of the TB 
mixing matrix, as is clear by comparing EqJTJJ to EqJTJ 

In the following discussion it is convenient use the expansion about TB mixing in- 
troduced in 1 19 1, 



r 1 1 

513 = 71' Sl2 = V3^ 1 + s ^ 523 = ~^ l + a ^ ^ l2 ' 

where the three real parameters r, s, a describe the deviations of the reactor, solar and 
atmospheric angles from their tri-bimaximal values. In terms of the deviation parameters 
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the la ranges for the mixing angles in EqsJHl El then translate into 
0.14 < r < 0.24, -0.05 < s < 0.02, -0.04 < a < 0.10, 



(13) 



with a central value of r = 0.2. To first order in r, s, a the lepton mixing matrix can be 
written as |19|, 
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Other related proposals to parametrize the lepton mixing matrix have been considered 
in [20]. 

In terms of the deviation parameters r, s, a, to leading order the TBR ansatz intro- 
duced in Eqs ]TTI]fTT1 is equivalent to having s = a = 0, corresponding to 



Sl3 



71' 512 "TT 523 "7T 

and a corresponding leading order mixing matrix, 
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which is equivalent to EqJTTJ with \U t 



e 3| = 0.14 corresponding to r = 0.2. By contrast, 
to leading order in r, the tri-maximal proposal [TBI H3, HE] corresponds to having s = 
and a = — (r/2) cos 5, which for r = 0.2 would imply, a = —0.1 cos<5. 

Having postulated the TBR form of the lepton mixing matrix, the next question is 
what neutrino mass matrix does this correspond to? In the flavour basis, in which the 
charged lepton mass matrix is diagonal and the TBR arises from the neutrino sector, 
the effective neutrino mass matrix {M^^) TBR is given in terms of the neutrino masses 
m 1 ,m 2 ,m 3 by, 

{M v eff ) TBR = C/ TBR diag(mi, m 2 , m 3 )U^ BK 

= mi$i$f + m 2 $ 2 $2 + m 3 $ 3 $3, (17) 
where we have written the mixing matrix in terms of three column vectors 

^TBR = ($l,$2,$3)- (18) 

It is clear that each column of the TBR mixing matrix in Eq{TU] can be written as a 
sum of a column of the TB mixing matrix in EqJT] plus a correction proportional to the 
reactor parameter r. Then, to first order in r, we find that the matrices in Eq{T7] can 
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be written as a sum of the TB matrices plus a correction proportional to the reactor 
parameter r, 
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$ 3 <&f = 5 | 1 1 | + ^ | 1 | . (19) 
From above we may write (M"ft) TBR as the symmetric matrix, 




(M: ff ) TBR = I . del, (20) 



where, 



2 1 

-mi + -m 2 , 

1 1 „> / 1 1 \ v/l 

— -mi + -m 2 - re I -mi + -m 2 ) + re [ -m 3 



1 « /I 1 



c = -mi--m 2 -re I -mi + -m 2 J + re ^-m, 

1 1 1 « f 1 1 

" = -mi + -m 2 + -m 3 + re -mi - -m 2 

6 2 \6 S 

1 1 1 iS f 1 1 

/ = -mi + -m 2 + -m 3 + re I -mi + -m 2 
111 

e = ~c m i - o m 2 + o m 3- (21) 
o 6 2 

In the limit that r = 0, (M^)™ reduces to the TB neutrino mass matrix (M"jj-) tb , 
and the relations b = —c and d = f and — e = a + 6 — d emerge as the characteristic 
signatures of the TB neutrino mass matrix in the flavour basis, in the convention for 
the TB matrix in EqJTl This implies that the origin of the reactor parameter r is due 
to a violation of the family symmetry that would lead to TB mixing. 

We now show how TBR mixing can arise at leading order from see-saw models based 
on sequential dominance (SD) [21J. To set the notation, recall that, in the type I see-saw 
mechanism, the starting point is a heavy right-handed Majorana neutrino mass matrix 
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Mrr and a Dirac neutrino mass matrix (in the left-right convention) Mp, with the light 
effective left-handed Majorana neutrino mass matrix M"ff given by the type I see-saw 
formula |22|, 



M v eff = M D M^ R M T D . (22) 
In a basis in which Mrr is diagonal, we may write, 

Mrr = diag(M A , M B , M c ) (23) 
and Md may be written in terms of three general column vectors A, B, C, 

M D = (A,B,C). (24) 
The see-saw formula then gives, 

AA T BB T CC T , . 

As noted in [10] {M^j) TB may be achieved if 

4 = -= I -1 V B = -= I 1 I . C = -= I 1 I. (26) 






which was referred to as form dominance (FD), since the constraints in Eqj26]lead to a 
form diagonalizable (M^j-) TB diagonalized by Utb (in this basis) with physical neutrino 
mass eigenvalues given by mi = a 2 /Ma, m 2 = b 2 /M B , m 3 = c 2 /M c . It is interesting 
to compare FD to Constrained Sequential Dominance (CSD) defined in In CSD 
a strong hierarchy \m-\\ |m 2 | < |m 3 | is assumed which enables ni\ to be effectively 
ignored (typically this is achieved by taking Ma to be very heavy leading to a very light 
mi) then CSD is defined by only assuming the second and third conditions in Eqj26] [llj. 
Thus CSD is seen to be just a special case of FD corresponding to a strong neutrino 
mass hierarchy. FD on the other hand is more general and allows any choice of neutrino 
masses including a mild hierarchy, an inverted hierarchy or a quasi-degenerate mass 
pattern. 

Clearly to achieve TBR mixing we must relax one or more of the conditions in 
Eqf26l There are many ways to do this, and here we just consider one simple example. 
Let us begin by considering a strong neutrino mass hierarchy \mi\ |m 2 | < \m^\ 
corresponding to Ma being very heavy leading to a very light mi as in CSD, in which 
the first condition in EqlSU] is irrelevant. As in CSD, we shall continue to assume the 
second condition in Eqj26]is accurately maintained, while we shall allow a small violation 
of the third condition in Eql26] parametrized by a small parameter e, 



B = ^\ 1 , C=^= 1 . (27) 
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We refer to this as Partially Constrained Sequential Dominance (PCSD), since one of 
the conditions of CSD is maintained, while the other one is violated by the parameter 
e. Note that the introduction of the parameter e also implies a violation of FD since the 
columns of the Dirac mass matrix B, C can no longer be identified with the columns of 
the MNS matrix, due to the non-orthogonality of B and C . In other words the R matrix, 
the general orthogonal matrix introduced by Casas and Ibarra [23], will differ from the 
unit matrix by an amount of order e, allowing non-zero leptogenesis as discussed in 

Assuming PCSD as in Eqj27l with a strong neutrino mass hierarchy \mi\ <C \ni2\ < 
1 777.3! , leads to a neutrino mass matrix, 



{K )PCSD 



BB T CC T 
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(29) 



To leading order in \m2\ 
(M v eff ) TBR in EqsflTIin] where we identify, 

m 1 = 0, m 2 = b 2 / Mb, m 3 = c 2 /M c , 



the mass matrix {M^j) PCSD in Eqs J28ll29l corresponds to 



re 



(30) 



Thus, the TBR form of mixing matrix in EqJT6lwill result, to leading order in |m2|/|m3|. 

It is straightforward to implement the above example of PCSD into realistic GUT 
models with non-Abelian family symmetry spontaneously broken by flavons which are 
based on the CSD mechanism [TTJ [T21 H3]- In such models the columns of the Dirac 
mass matrix in the diagonal right-handed neutrino mass basis are determined by flavon 
vacuum alignment, with the column B identified with a triplet flavon ^123 and the 
column C identified with a triplet flavon ^23 and it is quite easy to obtain a correction 
to the vacuum aligmment such that 



>123 




>23 




(31) 



in direct correspondence with Eqf271 For example, in such models based on the discrete 
family symmetry [13], the flavon vacuum expectation value (VEV) (0123) will preserve 
a Z2 subgroup of the original discrete family symmetry corresponding to an A4 generator 
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S |14j . while the flavon VEV (023) will violate this subgroup even in the limit that e — 0. 
It is therefore natural to assume some misalignment of (023) since, unlike ((^123) , it is 
not protected by any symmetry. 

To be concrete, in the A4 Pati-Salam model in [13] a radiative symmetry breaking 
mechanism is used, and in the first stage of symmetry breaking the alignment is achieved 
via an A4 invariant term but S'0(3)-breaking term A0^0*0^0* (summed over % = 1, 2, 3) 
which aligns the VEV of the field in the following possible directions, 

(|0|) oc (1,1,1) and/or <|0|) oc (1,0,0), (0,1,0), (0,0,1) (32) 

where only the magnitudes of the components are so far specified. What matters is 
the sign of the SO (3)-breaking term: if A > the "isotropic" option (|0|) oc (1,1,1) is 
picked up while the VEV is maximally "anisotropic" (i.e. with just one nonzero entry in 
(0)) if A < 0. Following this approach, in the first stage, triplet flavons 0123, 0i and 03 
are introduced and the above mechanism is used to align the magnitudes of the VEVs 
as 

<|0i23|> oc (1,1,1) and (10x1) oc (1,0,0), (|0 3 |) oc (0, 0, 1). (33) 

In the second stage, two further triplet flavons 023, 023 were introduced [13] and a poten- 
tial was assumed to lead to the CSD alignment. Here we show how a small modification 
of this potential can lead to alignment along the direction in EqJBTl We continue to as- 
sume that their VEVs are radiatively driven, but assume the slightly different alignment 
terms 

Am^S-^sl 2 + A1I0I.023I 2 + A1I0I.023I 2 + A23I023-023I 2 - (34) 

If A123 is positive, the VEV of 023 is driven to be orthogonal to (0123) while Ai > 
makes its first component vanish and thus (|023|) oc (0,1,1). In order to obtain the 
alignment of (023) in EqEH we now assume negative A23, which tends to align (023) 
along the direction of (023), and also assume that Ai is negative, which will tend to 
align (023) along the direction of (0i). The combined effect of these two terms will be to 
lead to an alignment of (023) of the form assumed in Eql3~Tl The message is that, since 
no symmetry is determining the alignment of (023), the choice of vacuum alignment 
for this flavon is not protected, and slightly different dynamical assumptions will lead 
to slightly different alignments. By contrast the alignment of (0123) is determined by 
simpler dynamics driven by the A4 symmetry. 

In such realistic GUT models there will also be additional corrections from TBR 
mixing due to charged lepton corrections, renormalization group (RG) running, and 
canonical normalization (CN) effects, as fully studied in |26j. In such a realistic frame- 
work, the TBR mixing matrix described here would correspond to the leading order 
neutrino mixing matrix at the GUT scale, with modified neutrino mixing sum rules as 
discussed in [26]. In such a GUT-flavour framework, one expects the charged lepton 
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corrections to the neutrino mixing angles to be less than of order 9f 2 /\^2 (where typi- 
cally 9\ 2 is a third of the Cabibbo angle) plus perhaps a further 1° from RG corrections. 
Thus such theoretical corrections cannot account for an observed reactor angle as large 
as 8°, starting from the hypothesis of exact TB neutrino mixing. 

In conclusion, we have proposed an extension of tri-bimaximal mixing to include 
a non-zero reactor angle #13 while maintaining the tri-bimaximal predictions for the 
atmospheric angle #23 — 45° and solar angle Q\ 2 = 35°. The TBR mixing proposal, which, 
to leading order in r, predicts the deviation parameters s = a = for all r, is distinct 
from the tri-maximal mixing proposal which, to the same approximation, predicts s = 
but a = — (r/2) cos 5. We have shown how such tri-bimaximal- reactor mixing can arise at 
leading order from the (type I) see-saw mechanism with partially constrained sequential 
dominance. We have in turn shown how partially constrained sequential dominance can 
be realized in GUT models with a non-Abelian discrete family symmetry, such as v4 4 , 
spontaneously broken by flavons with a particular vacuum alignment. In such models 
there will be theoretical corrections to the TBR mixing matrix due to charged lepton 
corrections, RG running and CN effects, but these are all expected to be subdominant 
compared to a large reactor angle # 13 in the region of 8°. Indeed, if RG/CN effects are 
very small and 9\ 2 is much smaller than the Cabibbo angle, then the tri-bimaximal- 
reactor predictions predictions for the atmospheric angle #23 = 45° and solar angle 
6 12 = 35° would be realized very precisely. If the present indication for 13 holds up, it 
will be very interesting to see how close the solar and atmospheric angles are to their 
tri-bimaximal values. 
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